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I. INTRODUCTION 

The understanding of Charge Dependence of strong in- 
teractions has been a crucial issue in Nuclear Physics (for 
reviews see e.g. [U-IIj])- In fact, the simplest place where 
this issue can be studied is the nucleon-nucleon interac- 
tion. As it is well known, isospin invariance is not an 
exact symmetry of strong interactions. As a consequence 
nuclear forces have a small but net charge-dependent 
component. By definition, charge independence means 
invariance under any rotation in isospin space. A vi- 
olation of this symmetry is referred to as charge de- 
pendence or charge independence breaking (CIB) and it 
means in particular that, in the isospin T = 1 state, the 
proton-proton (T3 = +1), neutron-proton (T3 = 0), or 
neutron- neutron (T3 = — 1) strong interactions are differ- 
ent. A particular case, known as charge symmetry break- 
ing (CSB), only considers the difference between proton- 
proton (pp) and neutron-neutron (nn) interactions. Fur- 
ther corrections are expected when, in addition, Coulomb 
forces are added to the proton-proton system (pp(c)). 

But, what is the scale of charge symmetry breaking?. 
Actually, the effects are important in the s-wave channel 
where an unnaturally large scattering length, due to a 
virtual state 1 in that partial wave, triggers a high short 
distance sensitivity. This, of course amplifies effects re- 
lated to variations in the short distance parameters pre- 
cisely in the region where the interaction and hence the 
charge symmetry breaking effects may be less reliable. 
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1 That means a pole in the second Riemann sheet in the negative 
energy axis. 



The current understanding is that CIB, and in partic- 
ular CSB, are due to a mass difference between the up 
and down quarks and electromagnetic interactions. On 
the hadronic level and in a One-Boson-Exchange (OBE) 
based picture (such as e.g. major causes of CIB and 
CSB are effects explicitly related to 

• Different proton and neutron masses. 

• Electromagnetic effects (mainly Coulomb interac- 
tion) . 

• Mass splitting of isovector mesons 7r and p and dif- 
ferent coupling constants. 

• Mass splitting between different A-isobar charge 
states. 

• Unknown short distance effects which are usually 
described by models. 

Traditionally it is believed that the difference between the 
charge and neutral pion mass in the One-Pion-Exchange 
(OPE) potential accounts for a big part of CIB while 
the difference between the masses of neutron and proton 
represents the most basic cause for CSB. Pion mass dif- 
ferences were shown to account for a 80% of the nn-pp 
scattering length difference |(|. The nucleon mass split- 
ting also generates a difference in the kinetic energies. 
Some recent OBE models only consider the differences 
coming from nucleon mass splitting and kinematical ef- 
fects However, these effects can only explain about 
a 15% of the empirical CSB. As a consequence some mod- 
els leave CSB unexplained Q while others simply intro- 
duce a term ad hoc to explain the remaining contribu- 
tion Q. In Ref. @ 27r-exchange contributions, irp di- 
agrams and other multi-meson exchanges including the 
A-isobar as intermediate states were considered to ex- 
plain the empirical CSB value accurately. In Ref. [l(| 
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27r-exchange contributions with A were found to be no- 
ticeable to explain the empirical CIB value being 3ir- 
and 47r-exchanges negligible. The difficulties arising in 
multi-meson exchange diagrams, in particular the energy 
dependence that they create, were avoided in the Bonn 
potential [B[ by introducing two effective scalar-isoscalar 
cr-mesons simulating 2-7T + irp exchanges. In the highly 
successful CD-Bonn potential CSB was included at 
the simplest one-boson-exchange diagrams with the same 
philosophy as its predecessor @- 

Many authors have also proposed the p — uj mixing 
as a key ingredient to understand CSB [l2|, (H)]. In 
Ref. [12j the p — uj mixing is identified as the major 
source of CSB while proton and nucleon mass differ- 
ences are identified to produce a minor effect. It should 
be noted that such a calculation is hampered by the 
fact that the g^NN coupling constant occurring in the 
CS part is about 40% larger than expected from SU(3) 
(g^NN — 3g P NN ~ 9) and also from the actual value 
taken for the CSB potential. Fixed-s Dispersion rela- 
tions [14] yielded g U NN = 5.7 ± 2.0 and Vector Meson 
Dominance uj — > e + e~ decays prefer g^NN ~ 10. The 
r\— 7T° has been shown to be of some relevance as well [lj|. 

The high sensitivity has been a major motivation 
to pursue experimental determinations of the neutron- 
neutron scattering length by indirect methods (for a re- 
view see e.g. Ref. [16j and references therein). A re- 
cent measurement of the nn scattering length using the 
pi-d capture reaction, yields (see also [l7| for a review 
) o-nn = —18. 69(4) fm when corrected for magnetic in- 
teractions. The recent CSB analysis of the reaction 
dd —> ir°a [l8[ uses also the large g^NN constant. 

The purpose of the present work is to approach the 
problem from a renormalization point of view. While we 
consider long distance physics to be known and describ- 
able by non-relativistic potentials we use a physical low 
energy parameter such as the scattering length to encode 
the unknown short distance physics. For the Charge In- 
dependent case it was shown how a natural SU(3) value 
of gtoNN could be confortably taken if a certain renor- 
malization condition was imposed (l9l [20| . As we will 
discuss in much detail this poses a problem of finiteness 
in physical observables when connecting different chan- 
nels such as np, nn and pp (strong or Coulomb). We 
propose a short distance renormalization condition fea- 
turing charge independence which guarantees finiteness 
although an ambiguity arises. However, a natural choice 
of the renormalization condition works quite well when 
compared to measured or recommended values. While 
the traditional point of view outlined above tried to com- 
pute the scattering lengths, the Effective Field Theory 
(EFT) approach assumes that these scattering lengths 
are completely unrelated (2l| - |26l j. We pursue here the 
possible connection between them from a new perspec- 
tive which actually is in-between, combining both points 
of view. We assume one scattering length to be known 
and exploit the concept of short distance insensitivity to 
determine all other scattering lengths and phase shifts 



from the requirement of finiteness of the scattering am- 
plitude. 

The paper is organized as follows. In Section [TT| we 
motivate the use of renormalization without assuming 
previous knowledge from the reader as a useful prepara- 
tory material for further developments. In Section Hill we 
extend the results to the interesting case of long range 
Coulomb interactions as they appear in pp scattering. 
CSB interactions are studied in all separate np,nn,pp and 
pp(c) cases in Section IIVI where insightfull universal re- 
lations are found out. In Section [V] we propose a short 
distance connection where all the channels are correlated 
from the requirement of finiteness. A further interest- 
ing application has to do with the determination in Ap- 
pendix [A] of the Gamow- Teller matrix element needed 
in the pp fusion process from the np scattering length. 
Finally, in Section IVTl we come to the conclusions. 

II. THE STANDARD VS. THE 
RENORMALIZATION APPROACH 

A. The OBE potential 

In this section we briefly review the main ideas behind 
renormalization in coordinate space for the OBE poten- 
tials (for a more detailed account see e.g. Ref. [19]) since 
they are focal in what follows. To provide a comprehen- 
sive perspective we compare it with the more traditional 
viewpoint of regulating the singular meson-exchange po- 
tentials by means of the introduction of short-distance 
form factors. The crucial distinction lies in the sensitiv- 
ity to short-distance details: from the renormalization 
point of view we expect complete insensitivity to these 
details. On the contrary, a regularization procedure only 
guarantees the finiteness of the results. For definiteness, 
let us analyze as an illustrative example the phenomeno- 
logically successful 1 S'o one boson exchange (OBE) po- 
tential [I, nil 

n 2 m 2 p- m -*r n 2 p —m a r 
y( r ) — y7TNN' rl TT e ticNN c 

1 ' ' 16ttM^ r A-k r 

, 9uNN e J P NN m P e 

4tt r MM% r ' 1 ' 

where g a NN is a scalar type coupling, g^NN a pseudo- 
scalar coupling, g U NN is a vector coupling and J p nn is 
a tensor derivative coupling (see || for notation). We 
neglect for simplicity nucleon mass effects and a tiny r\ 
contribution. We take m n = 138MeV, M N = 939MeV, 
to p = 770MeV, = 783MeV and g„ NN = 13.1 which 
seem firmly established. The OBE potential, Eq. (TTJ) 
corresponds to a long distance expansion of the poten- 
tial. On the other hand, NN scattering in the elas- 
tic region below pion production threshold involves CM 
momenta p < p max = 400MeV. Given the fact that 
l/m u = 0.25 fm <C l/p max — 0.5 fm we expect heavier 
mesons to be irrelevant, and p and uj themselves to be of 
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marginal important. This naive expectation is, however, 
not fulfilled in the traditional approach 0, [ll[ . 

In the following we will make the approximation m p = 
m u , specially when making fits of the coupling parame- 
ters to the So phases. Under the previous approximation 
is convenient to define 



9uNN 



9u>N N 



f, 



pNN m p 



2M% 



(2) 



in such a way that the combined u) — p potential reads 



oNN 



f 2 

J pNN""p 



TO 2 e -m p r 



9uNN e 



-in 



%nM% 



in 



(3) 



The previous simplification is useful since it avoids cor- 
relations between g^NN and J p nn in the 1 Sq channel. 
For SU(3) values of g uNN ~ 9 or VMD uj e+e" 
e.m. decays g^NN ~ 10.5 and typical p-tensor values 



fpNN ~ 14 - 18 one has g* 



ujNN 



0-7. 



B. Regular solution 

In the traditional approach [E|,[ll[ the problem is essen- 
tially handled by solving the reduced Schrodinger equa- 
tion, which for the s— wave case reads 



- <(r) + M N V(r)u k {r) = k 2 u k {r) . 



(4) 



with k = y/MffE the center-of-mass momentum, Mn 
the nucleon mass, and V(r) the boson-exchange potential 
of Eq. ([I]). The Schrodinger equation is a second order 
differential equation and it has two linearly independent 
solutions. The physical solution is usually determined by 
the regularity condition at the origin, i.e. 



«k(Q) = 0. 



(5) 



This boundary condition for the Schrodinger equation 
implicitly assumes that we are taking the potential seri- 
ously all the way down to the origin 2 . 

The asymptotic behaviour of the reduced wave func- 
tion for r>l/ m n is given by 



Ufc(r) -> 



sin(fcr -I- So(k)) 
sin<5o(fc) 



(6) 



where So(k) is the s-wave phase shift. For the poten- 
tial described by Eq. (jTJ , the phase shift is an analytic 
function of k with branch cuts located at k = dzim^/2, 
±im a /2, etc. This means in particular that for momenta 
below the first branch cut, \k\ < m^/2, we can expand 



the phase shift by means of the effective range expan- 
sion fH 



k cot<5o(fc) = — — + TT^o^ 2 
a 2 



v 2 k 4 



(7) 



where ao is the scattering length, r$ the effective range 
and i>2 the shape parameter. 

The effective range parameters can be related with the 
expansion of the wave function in terms of k 2 , i.e. 

Ufc(r) =u Q {r) + k 2 u 2 (r) + ... , (8) 

where uq and u 2 obey the following equations 

-u' \r) + M N V(r)u (r) - 0, (9) 
-u' 2 \r) + M N V(r)u 2 {r) - u Q (r) , (10) 

subjected to regular boundary conditions, uq(0) — 
""2(0) = 0, and asymptotically normalized to 

u (r) -> 1- r/a , (11) 
u 2 (r) -> (r 3 - 3a r 2 + 3a r r) /(6a ) , (12) 

for r 3> 1/™^. With this normalization, the effective 
range rg is computed from the standard formula 



ro 



2 / dr 
'0 



(1 - r/cto) - u (rY 



(13) 



In the traditional approach [TT| everything is obtained 
from the potential, which is assumed to be valid for 
< r < 00. In practice, strong form factors are in- 
cluded mimicking the finite nucleon size and reducing 
the short distance repulsion of the potential, but the 
regular boundary condition is always kept. One of the 
problems with this point of view has to do with the 
fact that the 1 5'o scattering length is unnaturally large 
ao = — 23.74(2)fm, while the effective range is natural, 
ro = 2.77(4)fm (approximately twice the pion Compton 
wave length, ~ 2/m n ). This has dramatic consequences 
regarding the short distance sensitivity, as we will show 
below. 

A fit to the np averaged data of Ref. [7j in the 1 Sq 
channel yields two possible solutions, see Table HI Thus, 
we have two good incompatible fits. A remarkable aspect 
is the fact that the vector meson coupling constant is 
accurately well determined. Actually, if we assume that 
we have fitted the potential, Eq. ([T]), to reproduce ao, a 
tiny change in the potential V — » V + AV has a dramatic 
effect on ao, since one obtains 



Aa 



a 2 Q M N 



AV(r)u (r) 2 dr : 



(14) 



2 Of course, in a more conventional setup strong form factors ac- 
counting for the finite nucleon size should be included. We will 
argue below that they play a marginal role in the discussion of 
CSB. 



a quadratic effect in the large ao- As a result, potential 
parameters must be fine tuned, as can be deduced from 
the previous fits. Thus, despite the undeniable success in 
fitting the data, this sensitivity to short distances looks 
counterintuitive. 
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BC 


r c (fm) 


m CT (MeV) 


g a NN 


duiNN 


x 2 /dof 


ao(fm) 


r (fm) 


Regular solution-I 





498.2(7) 


9.488(11) 


7.94(2) 


0.480 


-23.737 


2.678 


Regular solution-II 





550.72(4) 


13.87(13) 


20.10 (24) 


0.674 


-23.738 


2.679 


Renormalizing 





490(17) 


8.7(6) 


0(5) 


0.289 


input 


2.672 



TABLE I: Fits to the 1 S'o phase shift of the Nijmegen group 0] using the OBE potential with a charge dependent OPE part. 
We take m n a = 134.97 MeV, m n + = 139.57 MeV, g A = 1.29 and f„ = 92.4 MeV We neglect the CSB coming from the 

p-meson and take m p = m w = 770 MeV fitting m CT , QctNN and g^NN- We use the value a np = —23.74 MeV as an input when 
renormalizing. 
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FIG. 1: Zero energy wave function for the singlet np 1 So chan- 
nel as a function of distance (in fm) and for the different sce- 
narios with large and small ui— couplings. This wave function 
goes asymptotically to Uo(r) — s> 1 — r/ao with cto — — 23.74fm 
the scattering length in this channel. The zero at about 
r = 0.5fm signals the existence of a spurious bound state. 



However it is worth mentioning that the different sce- 
narios correspond to selecting a potential possessing spu- 
rious bound states or not. The spurious bound state 
problem has been discussed in Ref. 19] at length; the 
number of inner zeroes of the zero energy wave function 
provides the number of bound states. For illustration we 
represent the zero energy wave in Fig [TJ In the regular 
case, the OBE potential with a big g^NN ^ s f ree °f bound 
states. However if a small S^ntv i s chosen, then one has 
to deal with a bound state which does not exist and it is 
hence spurious. 



C. Irregular solutions 

The previous results are in conflict with the intuitive 
expectation of insensitivity of low energy physical observ- 
ables with respect to the specific details of the potential 
in the short distance region. Otherwise, where should 
one stop ?. This is the basis of the renormalization view- 
point. The way to proceed is to impose renormalization 
conditions which eliminate the short range sensitivity at 
the expense of treating low energy parameters as inde- 



pendent variables from the potential. An example of a 
renormalization condition (RC) is to fix the scattering 
length, with the consequence of avoiding the fine tuning 
problem summarized by Eq. (fT4"|) . In other words, we 
trade the explicit dependence of the results on the short 
range parameters of the potential for low energy observ- 
ables. The values of the later are usually well-known by 
other means. 

In principle there are several ways in which one can 
impose renormalization conditions, one popular example 
being counterterms. They correspond to the coupling 
constants of a short distance contact potential which is 
expanded in terms of 5 functions and its derivatives 



V c {$) = Co 5(x) + C 2 {V 2 , 8(x)} + 



(15) 



where the dots represent terms involving higher deriva- 
tives of the 6 function. This potential is added to the 
usual finite range potential V , and then the correspond- 
ing Schrodinger equation for Vc + V is solved. The re- 
sulting potential is strongly singular and needs to be reg- 
ularized by introducing a cut-off r c , a length scale which 
is used to smear the S functions. The different coupling 
constants Co(r c ), C2(r c ), etc, are set to reproduce the de- 
sired low-energy observables. A nice presentation of the 
previous method is given in Ref. [29]. The disadvantage 
is that the procedure of using a potential to renormalize 
quickly runs into problems when one tries to decrease the 
size of the cut-off. For example, it may be impossible to 
reproduce certain physical observables, specifically the ef- 
fective range, if the short distance cut-off is too small |30| . 
unless one accepts complex values for the counterterms 
Co(r c ) and C2{r c ) which violate either causality or off- 
shell unitarity (see Ref. [3l[ for a detailed discussion). 
Related positivity conditions are discussed in Ref. [l9j . 

Here we use a more indirect method to renormalize 
which is able to avoid some of the previously mentioned 
complications. The idea is to substitute the regularity 
condition of the Schrodinger equation, Ufc(0) = 0, by an 
arbitrary boundary condition at the origin 



ifc(0) = 



<(Q) 

u fc (0) ' 



(16) 



where we have used the log-derivative of the wave func- 
tion, instead of an independent condition for Ufc(0) and 
tti(0), as that will only affect the precise normalization 
of the wave function, which can be later determined by 
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other means. The regularity condition u k (0) = cor- 
responds to taking the limit L k (0) — > oo (as u' k (0) is a 
constant), but by changing the precise value and energy 
dependence L k (0), the values of low energy observables 
can be fixed. 

The previous procedure would in principle involve a 
fitting procedure, which can be avoided by taking into 
account the expansion in powers of k of the wave func- 
tion, Eq. ([8]). For example, if we want to fix the scat- 
tering length, we will solve the corresponding equation 
for the zero-energy wave function Uo(r), with the asymp- 
totic (r — > oo) boundary condition of reproducing the 
scattering length 



u' \r)+M N V(r)u (r) = 0, 



uo(r) 



1 



r 

a ' 



(17) 
(18) 



but, instead of solving the previous equation from r = 
to r — » oo, we solve it downwards from infinity to the 
origin. Then, we assume that U2{r), it 4 (r), etc, are 
subjected to regular boundary conditions at the origin, 



u 2 (0) = u 4 (0) 
to take 



and 1*2(0) = u' 4 (§) = 0, which means 



ifc(0) 



«o(0) 
uo(0) 



(19) 



By so doing we achieve some insensitivity at short dis- 
tances, as we will show later 3 . Fixing more scattering 
parameters is straightforward: one solves downwards the 
corresponding equations for uq{t), 1*2 (r), . . . , U2 n (r) with 
the asymptotic conditions of reproducing ao, To, . . . , v n , 
and assumes trivial boundary conditions for U2 n +2(r), 
U2n+A(f)t e tc, resulting the following logarithmic bound- 
ary condition 



ifc(0) 



u (0) + fcV 2 (0) 



k 2n u' 2n (0) 



uo(0) + k 2 u 2 {Q) + ■■■ + fc 2n u 2 „(0) 



(20) 



In practical computations it is convenient to introduce a 
short distance cut-off, r c , and then take the limit r c —> 0. 

A further simplification can be made if we take into 
account that the OBE potential of Eq. ([I]) is local and 
energy independent. This means in particular that dif- 
ferent energy states are orthogonal, 



u k (r)u p (r)dr = , 



(21) 



for k 7^ p, which requires an energy independent bound- 
ary condition at the origin, as a consequence of the next 
equality 



u k {r)u p (r)dr = u' k u p - u p u' k 



(22) 



3 In dispersion theory these renormalization condition resembles 
the customary subtractions. In our case the form of the subtrac- 
tion is a bit more involved as discussed below. 



which means that the orthogonality condition Eq. 
can be re-expressed as 

<(0) _ "p(0) ; 
Ufc(0) Up(0) ' 



(ED 



(23) 



or, equivalently, L k (0) = L p (0), implying an energy in- 
dependent boundary condition. 

Here we will only consider the case in which orthogo- 
nality is preserved. The restriction is that orthogonality 
implies that we can only fix one scattering parameter, 
namely the scattering length. Therefore we will integrate 
downwards the zero energy state, Uo(r) — > 1 — r/ao, up 
to the cut-off r c . At this point, we can make use of the 
superposition principle in order to take real advantage 
of the boundary condition method. For the zero-energy 
solution, the superposition principle can be used to write 
the wave function as a linear combination of two inde- 
pendent zero-energy solutions 



u (r) = ui(r) - 



u r (r) 
a 



(24) 



where u\(r) and u r (r) are solutions of the zero-energy 
Schrodinger equation, Eq. ([9]), with the asymptotic 
boundary conditions 



ui(r) 
u r {r) 



(25) 
(26) 



at large distances. The previous expression for uq as a 
linear combination, Eq. (|24[) . can be introduced in the 
integral representation of the effective range ro , Eq. (jTJ 
yielding the following correlation between rg and ao 



r = A -\ h — , 

where Ao, Bq and Co are given by 

dr(l-uf), 

X) 

dr(r — uiu r ) , 
dr(r 2 - u 2 r ) . 





poo 










Bo 


= - 4 / 




Jo 




poo 


C 


= 2 




Jo 



(27) 

(28) 
(29) 
(30) 



The interesting feature is that the dependence of the ef- 
fective range with respect to short-range parameters of 
the potential is greatly diminished. 

The short distance sensitivity can be vividly seen in 
Fig. [21 where the regular (parabola like curve) as well as 
the renormalizcd (flat curve) effective range for the OBE 
potential are shown as a function of g^NN- For simplicity 
only the solution with the small <7^,jvat (Regular solution 
I) is represented. 

The finite-energy solutions and the phase shifts can be 
obtained from the orthogonality condition, which implies 



U 'k( r c) 

rT-^o u k (r c ) 



lim 



lim 



u o( r c) 

u Q (r c ) 



(31) 



6 



a 
c 




the following formulas 



g*coNN 



FIG. 2: Dependence of the effective range with respect to 
9u)NN i n the regular case with a small coupling constant and 
in the renormalized one. 



providing the initial boundary conditions for the finite- 
energy Schrodinger equation, Eq. (|4]). We normalize the 
scattering wave function as follows 



u k (r) 



sin(£;r + So) 
sin 60 



(32) 



Again, if we use the superposition principle, Uk can be 
written as 



Ufc(r) = k cot S u k ,s(r) +Ufc, c (r), 



(33) 



with u kc and u^ s solutions of Eq. Q asymptotically 
normalized as 



Uk.c — > cos(fcr) . 
sin(fcr) 



Uk,. 



(34) 
(35) 



These two wave functions have the property of going to 
u-k,s ~~ ^ u r and u kc — ► U\ for k — > 0. Then, the orthogo- 
nality constraint Eq. (|3"Tj) reads 



k cot 8 u' k s {r c ) + u' k c {r c ) _ g ?4(r c ) - <(r c ) 
kcot5ou k ,s(r c ) + u k . c (r c ) a ui(r c ) - u r (r c ) 



, (36) 



A ( U\ 


= lim (ui(r c )u' k (r c ) 


- u x 


B(k) 


= lim (u k , c (r c )u' r (r c ) 


— u r 


C(k) 


= lim (u[(r c )u k , s (r c ) 


- Ul 


V(k) 


= lim (u r (r c )u' k (r c ) 





(38) 

The previous expressions fix the arbitrary normalization 
of Eq. The obvious conditions .4(0) = D(0) = and 
S(0) = C(0) = 1 are satisfied. Expanding the expression 
for small k one gets that v k is a polynomial in l/oeo of 
degree k + 1 . 

Finally, we can use the previous procedure to fit the 
np averaged data of Ref. [7| in the 1 5o channel (once 
we have fixed the scattering length to its experimental 
value) , yielding the values in table fl] We can see the 
large uncertainty on the value of g^NN' which shows that 
there is a greater insensitivity to shorter distances after 
renormalizing. This agrees with the previous remarks 
on the sensitivity of the effective range on g^NN' illus- 
trated in Fig. [2j Let us note further that as discussed 
in Ref. the renormalization scenario also has also a 
spurious bound state as in the small <?^,jvAr regular solu- 
tion case (see FigQJ. The current discussion would be 
modified by the inclusion of form factors which incor- 
porate the finite nucleon size. However, because of the 
short distance insensitivity form factors turn out to play 
a marginal role (l9| after renormalization . 



D. Review of the Renormalization Process 

The renormalization procedure proposed in this section 
can be summarized as follows 

• For a given scattering length ao, integrate in the 
zero energy wave function uo(r) with Eq. down 
to the cut-off radius r c . This is the renormalization 
condition. 

• Implement self-adjointness at the cut-off radius 
through the boundary condition 



u' k (r c )u (r c ) - u' (r c )u k (r c ) = , 



(39) 



where we have dropped the lim rc _j.o for shortening the 
notation. Note that the dependence of the phase-shift 
on the scattering length is explicit: k cot So is a bilinear 
mapping of cvq, 



k cot So = 



a A(k) + B(k) 
a C{k)+V{k) 



(37) 



where the functions A, B, C and V are even functions of 
k which depend solely on the potential and are given by 



• Integrate out the finite energy wave function u k (r) 
with Eq. ((3]) and determine the phase shift So(p) 
from Eq. ©. 

• Remove the cut-off (take the limit r c — > 0) to assure 
model (regulator) independence. 

This allows to compute Sq (and hence ro, «2 ) from (i) 
the potential V(r) and (ii) the scattering length ao as 
independent information. Note that this is equivalent 
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to consider, in addition to the regular solution, the ir- 
regular one. In momentum space this can be shown to 
be equivalent to introduce one counterterm in the cut- 
off Lippmann-Schwinger equation (see Ref. [3l| for a de- 
tailed discussion). Both Eq. ([27]) and Eq. ([37]) highlights 
this de-correlation between the potential and the scatter- 
ing length. Contrary to common wisdom, but according 
to our intuitive expectations, no strong short range re- 
pulsion is essential. The moral is that building a from 
the potential is equivalent to absolute knowledge at short 
distances, as in the 1 Sq channel a strong fine tuning is at 
work. This example illustrates our point that the renor- 
malization viewpoint tells us to what extent short dis- 
tance physics may be less well determined than the tra- 
ditional approach assumes. This opens up a new perspec- 
tive (see |l9() to the phenomenology of OBE potentials 
in cases where the strong w-repulsion has proven to be 
crucial at low energies. 



III. RENORMALIZATION WITH COULOMB 
INTERACTIONS 

In this section we generalize the previously discussed 
renormalization approach to the case of proton-proton 
scattering, where the infinite range of the Coulomb in- 
teraction will pose some problems. The corresponding 
s— wave reduced Schrodinger equation is 



CI' 1 



M, 



,c 



k 2 u c k . 



(40) 



where k is the center-of-mass momentum, M p the pro- 
ton mass, V pp (r) the strong proton-proton potential, and 
a ~ 1/137 is the fine structure constant. Actually, the 
current discussion is tightly linked to the corresponding 
one for the two potential formula presented by two of us 
recently [321 ], 



(see for example [33[). At short distances these solutions 
behave as 



u £rcg( r ) 



a,B 3 a 2 



(44) 



J 0,irr 



V) -+ 1 



2r 



(IB 



log 



B 

2r 



2 lE -\ 



2r 

(IB 



2 r 



1 , 2r 5 
- log — + 7B - - 

I CLB 4 



+ 0(r 3 ), 



(45) 



where = 0.57722 is the Euler-Mascheroni constant. 
The previous means in particular that Vq reg is the short 
distance regular solution and Vq irr the short distance 
irregular solution. In principle, in the absence of any 
strong potential, Vq r would be the zero energy solution 
for the repulsive Coulomb potential. The presence of the 
strong interaction between the protons means that the 
zero-energy asymptotic solution for r — ¥ oo will be in 
general a linear combination of i>oJ ra g and v§ irr . 

For the proton-proton system the Coulomb scattering 
length is related with the asymptotic behaviour of the 
zero energy wave function at large enough distances 



u 0,reg 



('■) 



ao,c 



(46) 



where i>oJ reg and Vq irr are the previously defined regular 
and irregular zero energy wave functions, and ao,c is the 
s— wave Coulomb scattering length. If the Coulomb in- 
teraction is switched off by taking as — > oo, the previous 
wave function reduces to the corresponding one for finite 
range forces, v(r) = 1 — r/ao. 



B. Effective range 



A. Coulomb scattering at zero energy 



The Coulomb effective range is given by the following 
formula 



The longest range piece of the proton-proton inter- 
action is the Coulomb repulsion between the protons. 
Ignoring any strong effects, zero energy proton-proton 
scattering in s— waves can be described by the reduced 
Schrodinger equation 



a sr 



(41) 



where as is the proton Bohr radius, which is defined as 
as = 2/MpCt = 56.62 fm. The previous equation has the 
following two linearly independent solutions 



y 0,irr 



(r) = 2v^Xi(2VS), 



(42) 
(43) 



where x = 2r/as and Ki(x) and I\(x) are modified 
Bessel functions of the first and second kind respectively 



r , c = 2 / *[, c (r) 2 - Uo c (r) 2 ] , (47) 
Jo 

where Vq (r) is the Coulomb zero energy solution given 
by Eq. (|46]) , and Uq (r) is the full zero energy solution to 
the Schrodinger equation 



uf+( M p V pp (r) + ^ z ) u° (r) = , (48) 



2 

a B r 

subjected to the asymptotic boundary condition 

"oW-ntfW, (49) 

for r — > oo. By making use of the superposition principle, 
the solution Uq can be decomposed as 



Uoi r ) = u 0,h-r( r ) ~ °' reS ^ ^ 



(50) 
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where u^ hr and w^ rcg are solutions of the zero energy 
Schrodinger equation, Eq. (|48p , behaving asymptotically 
(r — > oo) as 



^O.irr 



;( r ) ~> ^regW> 

0) -> ua irr (r). 



(51) 
(52) 



The subscripts rcg and ; rr do not refer to the regularity 
of the solutions at the origin, but with the long range 
behaviour of the full solutions. 

By plugging the decomposition of the full and purely 
Coulomb wave functions, Eqs. (|4"!)f and (dHJ), into the 
integral representation of the Coulomb effective range, 
Eq. (|47|). we obtain the following correlation between the 
Coulomb scattering length and effective range 



ro,c 



A° 



-°o 



eg 



l o,c 



(53) 



which is a direct generalization of Eq. (|27[) for the non- 
Coulomb case. The Coulomb correlation functions Aq , 
Bq and Cq are given by the integral expressions below 



A° = 2 / dr{ V C(rf-uC{rf) 



(54) 



5 ° = -4 / ^(«0,irrW<regW-<irrW«0,regW), 
JO 

(55) 

C C = 2 / dr« rcg (r) 2 -^ rcg (r) 2 ). (56) 
Jo 



C. Coulomb scattering at finite energy and 
Coulomb effective range expansion 

The definition of the phase shifts in the presence of the 
Coulomb potential is related with the behaviour of the 
wave function at long distances, which is given by 



?4 (r) -> cot 6% (k)F (r), p) + G (»/, p) 



(57) 



where %{k) is the Coulomb- modified proton-proton 
phase shift, k the center-of-mass momentum and Fq(t], p) 
and Go(?7, p), with rj = 1/fca.B and p = kr, are the 
s— wave Coulomb wave functions [33]. The wave func- 
tion is the solution to the reduced Schrodinger equation 
Eq. (|4T)]) . The Fo(rj,p) and Go(r],p) wave functions be- 
have asymptotically (r — > oo) as 

F -> sin (At - r? log(2fcr) + cr ) , (58) 
Go -> cos(fcr - ?7log(2A;r) + cr ) , (59) 

with (To the Coulomb phase shift, which is defined as 



r(i-wj) 



(60) 



The phase shift in presence of the infinite-ranged 
Coulomb force does not obey the usual effective range 



expansion, valid for short-ranged potentials, but a 
Coulomb-modified effective range expansion, given by 



k cot# C\n) + — h{rj) 



with G{rj) and /i(r?) defined as 



1 



1 



oo 

£ v n ,ck 2n 



ro,ck 



c z (v) 



e 2ir V _ I 



(61) 



(62) 



1 



^— ' 7i(n 2 + n 2 ) 

n— 1 x ' 



\ogi]-j E . (63) 



For obtaining the Coulomb extension of Eq. (f3"T|) we 
use again the superposition principle to write lijj? in the 
following way 



Cfa) ?4 (r) = [k cot (5 c C 2 ( V ) + — h(v) < reg (r) 

- ^irrW, (64) 

where u^ reg and u^ ilT are solutions of Eq. (I^ . which 
obey the asymptotic boundary conditions 



U fc,reg( r ) -> 



(65) 



U £ lrr (r) H- -C(»/)Go(»7,p) + ^)Fo(»j,p), 

C(?7) 

(66) 

for r — > oo. These two solutions have been normalized 
in such a way that for k — > they coincide with the 
previously defined zero-energy wave functions UQ YCg and 
ti£ irr , see Eqs. flST) and ([52]) . 

Once these definitions have been made, it is straight- 
forward to obtain the correlation 



k cot 6$C 2 {rj) + — fcfa) 



«o.c-4 c (fc) + g c (fc) 
a , c C c (k)+V c (k) ' 



(67) 



where -4 c (fc), B°(k), C c (k) and X> c (fc) are defined as 



A c (k) = lim 
r ->-Q 



= lim 

r o ->0 



C c (fc) = lim 



X> G (fc) = lim 

TV— >0 



u C^irr( r c) u fe,iii ( r c) 
U 0^irr ( r c) M fcarr( r c) 
M arcg( 7 'c)Ufe i i I . r '(r c ) 

M fe!irr( r c) M ?reg ( r c) 



u 0,irr( r c)w fc reg (r c ) 



U 0,irr ( r c) w fe,reg( r c) 
U 0,rcg( r c) u k,rcg ( r c) 
U 0,reg ( r c) u fc,reg( r e) 



(68) 



(69) 



(70) 



(71) 
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IV. CHARGE SYMMETRY BREAKING 

In the previous sections we have shown how the renor- 
malization of the two nucleon system can be carried 
out. This procedure allows to determine the 1 So phase 
shifts for np, nn, pp and pp{c) from their corresponding 
scattering lengths a np , a nn , a pp and a pp respectively. 
The previous computation can be compared with the ex- 
perimental values for these quantities in order to test the 
renormalization procedure. 

Admitted values of the scattering lengths are 



^0,pp 
&0,nn 



-7.8149(29) fm, a , pp 
= -18.8(3) fm , a ,n P 



-17.3(4) fm, 
-23.77(9) fm, 



(72) 



giving AaciB = 5.7 fm and AacsB = 1-5 fm. For the 
effective range we have 



„c 



2.769(14) fm, 
= 2.75(11) fm , 



' pp 



2.85(4) fm, 
: 2.75(5) fm, 



(73) 



with Atq^cib = 0.05 fm and Ar^css = 0.1 fm. As can 
be seen, the CIB/CSB breaking is much larger for the 
scattering length than the effective range. Part of that 
is explained by the unnaturally large value of the NN 
scattering length. 

To take into account the various physical effects which 
generate charge symmetry breaking, we consider the 
neutron-proton mass difference and the OPE potentials 



^ nn V / 



with 



-M v fl 



-M np fl 



2 



r 

-m^ n r 



171^ + 



(74) 



: + = 139.57 MeV. M np is 



„ = 134.97 MeV and m 
twice the np reduced mass, 2fi np — 2M p M n / (M p + M n ) 
Therefore, for the OBE AW potential we have 



V np (r) 
V nn {r) 
V PP (r) 



V la (r) 



(75) 



Clearly, the potentials in the different channels are not 
very different from one to another quantitatively. Actu- 
ally, the a and uj exchange contributions coincide identi- 
cally. On the other hand the 7r and p take into account 
the different charged mesons which are exchanged. Obvi- 
ously, one expects the symmetry breaking effects coming 
from 7r exchange to be more important than p exchange. 



Theoretical computations seem to support the previous, 
giving Ao:cib,tt = 3.24 fin and AaciB.p = —0.29 fm. see 
Ref. [l(|. As a consequence p mass differences are negli- 
gible. 

The long distance correlation between the scattering 
length and effective range looks as 



>0,np 
^0,pp 

r C 
0,pp 



A n 



A n 



B,, 



«0,: 

B. 



+ 



On 



0,np 



PP 



-Ann 



Brill 



0,np 



a 



rip 



CtO,nn 

B c 

A c , D VP 

pp ' 



0,np 



pp 



OtO,C,pp a 0,C,pp 



while the phase shifts are given by 



k cot So,np 
k COt 5o,nn 

kcotd 0ynp 
C 2 (n)fc cot 5% pp 



ao,npA n p(k) + B np (k) 

Ceo,npC n p(k) + T> np {k) 

A„ n {k) + B nn (k) 

a0,nrAm(fc) + ^>nn[k) 

ao,npA n p(k) + B pp (k) 
ao,ppCpp(k) + V pp {k) 

— h(n) 
a B 

a o, PP A pp (k) + B pp (k) 



0,pp pp 



Cg(k)+Vg(k) 



(76) 
(77) 
(78) 
(79) 

(80) 
(81) 
(82) 

(83) 



We remind that the scattering lengths are independent of 
the potentials. 

In Fig. [3] we show the universal functions A, B C and 
T> for the four cases considered. As can be seen, for nn, 
np and pp they coincide even though the potentials are 
different. This means in particular that most of the CIB 
and CSB effects for p < 400MeV come solely from the 
difference in the scattering length (there are no genuine 
sizeable effective range effects). It is also interesting to 
see that the Coulomb corrections to the pp(c) universal 
functions differ increasingly for higher energies. 



V. THE SHORT DISTANCE CONNECTION 

As it is well known, at large energies the nn, np and 
pp phase-shifts start to resemble each other, which means 
that charge invariance is respected for large enough mo- 
menta. Most of the charge invariance and charge sym- 
metry breaking effects only affect the low energy be- 
haviour, specifically the scattering lengths, where one 
finds Aacm = 5.7(3) fm and AacsB = 1.5(5) fm. When 
one considers the effective range, the symmetry breaking 
effects are already ten times smaller than in the scatter- 
ing length case, being of the order of the tenth of a fermi. 
The problem is how to explain these differences. 
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FIG. 3: The universal functions A, B C and T> defined by Eqs. (|68I71[) in appropriate length units as a function of the CM 
momentum p (in MeV) for the four 1 So channels np,pp,nn,pp(c) . These functions depend on the potentials V np (r), V nn (r), 
V pp (r), V p c p (r) only but are independent of the scattering lengths. 



In the traditional approach all the CIB and CSB ef- 
fects are explained via the OBE potential, Eq. (TTJ). The 
Schrodinger equation is integrated from the origin to in- 
finity with regular boundary conditions and all the differ- 
ence between scattering observables must come from the 
potential. In the renormalization approach things get 
more involved: there are explicit contributions coming 
from short distance operators which are used to weaken 
the short distance sensitivity. The problem is how to 
implement either charge independence or its breaking 
within this approach in a regulator independent way. If 
we assume that at lowest order all the charge indepen- 
dence breaking comes from the finite range potential, one 
is tempted to identify short distance charge independence 
with identical logarithmic boundary conditions. For ex- 
ample, if we relate the nn and np problems with 

— r~\ ~ — ' ( ' 

we will find that this relation produces log-divergent re- 
sults for the OBE potential in the limit r c —> 0. Another 
option is to regulate with a short distance delta potential 

Vc(r;r c ) = ^^5(r-r c ), (85) 



which corresponds to a specific regularization of the 5 
function potential, and assume that charge independence 
at short distance is equivalent to Co, nn (r c ) = Co^ np (r c ) = 
Co,pp( r c) = Co(r c ). This choice leads to the following 
logarithmic boundary condition between nn and np 

1 f u' nn {r c ) l\_ 1 f u' np {r c ) 1 \ 

where M n is the neutron mass and M np is twice the re- 
duced np mass. The counterterm conditions also runs 
into the same cut-off dependence problems than the log- 
arithmic boundary condition. This means in particular 
that the two previous proposals are regulator dependent, 
and hence model dependent, and pose a serious problem 
on what is meant by charge independence of short dis- 
tance operators. We will show that by using the hypoth- 
esis of charge independence at short distances together 
with finiteness, a relation between them can be estab- 
lished which works rather satisfactorily. 

At short distances all the pp (strong/Coulomb), np and 
nn potentials have an attractive Coulomb like behaviour 

V NN (r) ► — — , (87) 

r-s-0 Hr 

where NN either refers to pp (strong), pp (Coulomb), 
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np or nn, and [Xnn and Vnn are the corresponding re- 
duced mass and potential. The constant R depends on 
the problem; for the OBE potential of Eq.([T]) with the ad- 
ditional simplification of taking — m p and defining 
9luNN ' we get tne scales 

— — = M np (fl NN + gl NN - g^NN 2 ) > ( 88 ) 
— — = M n (fl NN + gl NN - g^NN 2 ) > ( 89 ) 
-L = M p (fl NN + g 2 aNN - g* uNN 2 ) , (90) 

pC = M P (f 2 NN + glNN - glNN 2 ~ a ) > (91) 



with M np twice the reduced np mass, M n 



' l^np- 



As 



a consequence of the short distance Coulomb singularity, 
the wave function at short distances approximately be- 
haves as linear combinations of attractive Coulomb wave 
functions 

R 

uk,NN{r) -> a — y/xJi(2y/x) + b2y/xYi(2'</x) 



+ 0(mr,mR, 



(92) 



where now the Bessel functions J\ and Y\ are used (in- 
stead of 1\ and K\). The constants a and b determine 
the correct linear combination, R can either be i? nn , R n p, 
and i? pp (strong/Coulomb), x — 2r/R and m generically 
denotes the mass of any of the exchanged bosons. The 
expected mR contributions will only shift the irregular 
solutions by a constant. 

The previous behaviour can be quite problematic as 
we can see if we consider the log-derivative of the wave 
function at small enough cut-off radii, which behaves as 



R 



k.NN 



(To) 



Uk,NN(r c ) 



-27- 



-RX 

4 



log — 

B R 



(93) 



where A = ag/bo and the dots refer to higher order terms, 
like mr c or k 2 r^ corrections. With this behaviour, we 
can see that naively identifying the log-derivative at the 
cut-off radius in order to obtain correlations between ob- 
servables of the different two nucleon systems will yield 
divergent results. For example, relating np and nn 



k,np 



(rc) 



^fc,nn(^c) ^fc,np v c) 

generates the singularity 



log [ — — 



i 



Rn 



■ loe 



Rn 



(94) 



(95) 



This singularity is indeed mild, as it can only be seen 
at very short distances (depending on how small is the 
difference between 1/R nn — 1/-R. 
will ruin our results. 



np)i 



but sooner or later 



Under these circumstances there is a quantity that can 
be constructed from the log-derivative at short distance 
that is finite in the r c — > limit. This quantity is the 
following 



S = R 



u'(r c ) 
u(r c ) 



(96) 



which is cut-off and energy independent. This suggests 
that different scattering problems, having different short 
distances constants but the same logarithmic scale de- 
pendence, can be connected in such a way that the scale 
dependence is eliminated. This is done by equating the 
corresponding 5's 



S\ — 6>2 



(97) 



where 1 and 2 refer to two different NN — nn, np, pp, 
pp(c) cases. We can give here two examples of the ad- 
equacy of the short distance connection. The first one 
is to obtain the strong pp scattering length from the ex- 
perimental Coulomb one, a$ pp — — 7.8149fm yielding 
Q^pp = -18.46 fm, a not unreasonable results (to be com- 
pared with the extraction ao :PP — — 17.3 fm, see Ref. 0], 
where the error comes from model-dependence). The 
CD-Bonn potential gives a value of ao tPP = — 17.46 fm. 
The extracted effective ranges are Tq pp = 2.735 fm and 
; 2.789 fm. As a second example, by taking the 



r , P: 

np scattering length as input, ao, 



-23.74 fm, we 



can obtain all the NN scattering lengths, giving «o,nn = 
-19.626 fm, a . pp = -17.806 fm and a% pp = -7.706 fm 
for the scattering lengths and ro, n p = 2.672 fm, ro, n n 
■2.771 fm, r ,p P = 2.802 fm and r% pp = 2.747 fm for the ef- 
fective ranges. A remarkable aspect of the previous com- 
putation is that one obtains AaciB = 5.024 fm, Arcm = 
0.115 fm, Aa CSB = 1.82 fm and Ar CS B = 0.031 fm which 
agree within error estimations with the expected values 
for these two quantities [2| . In Table [TT] we summarize 
the results obtained with the short distance connection 
(renormalized) and the one obtained integrating upward 
with a regular boundary condition (regular). We can 
see that in the case of a big g^^N l ne regular solution 
does a poor job in calculation the low energy parameters 
(LEP) in other channels. The CD-Bonn potential fu\ 
corresponds with this scenario, i.e., a big SU(3) breaking 
coupling constant but with any spurious bound state. 
Looking at this table one can understand why in this 
model a different mass for a ficticious er-meson is used in 
each NN channel. The strong fine-tuning that appears 
in this situation hinders the relation between different 
NN problems. 

A further interesting example of the adequacy of the 
short distane connection is illustrated in Appendix [A] 
where the Gamow- Teller matrix element appearing in the 
proton-proton fusion process is analyzed. 

Note, however, that the previous is not the only possi- 
ble covariant short distance connection, as we could have 
defined 



S' = R 



u'(r c ) 
u(r c ) 



(98) 
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FIG. 4: The relations between the predicted scattering lengths in the 1 So channel for nn,pp,pp(c) as a function of cto y n P when 
the successive n, n + a and n + a + p + uj contributions are included. We plot inverse scattering lengths. Note the small scale. 



with A some arbitrary constant, which depends on the 
specific NN problem which is being considered. A nat- 
ural choice is to take A of order unity, which does not 
make much difference between different choices of S due 
to the weak logarithmic behaviour. It must be stressed 
though that the results are not unique: arbitrary A's 
can be introduced to better connect the different two 
nucleon systems. As the hypothesis of charge depen- 
dence of short distance operator cannot be implemented 
in a completely model independent way, we will chose to 
take X nn = X np = X pp at first order. We have already 
seen that this simple condition generates quite accurate 
results, meaning that corrections due to AA are indeed 
small, and can be effectively considered as higher order 
effects, confirming thus our expectations. As an example 
of what values of A A to expect, if we try to correlate 
the strong and Coulomb scattering lengths, we will get 
X pp - Ag, = 0.0321 - 0.0471, where the range given is a 
consequence of the uncertainty in a pp = —17.3(4) fm. 

To clarify the implications of the short distance con- 
nection, let us consider two different problems 1 and 2, 
which have the associated Coulomb length scales i?i and 
R 2 . In other words, we have the differential equations 

-u'( + 2 t i 1 V 1 (r)u 1 (r) = k 2 Ul (r) , (99) 
-u'^ + 2^ 2 V 2 {r)u 2 {r) = k 2 u 2 {r) , (100) 

where the reduced potentials behave as 1/r at short dis- 
tances 



2fi 2 V 2 (r) 



1 



Rir 
1 



(101) 
(102) 



These two problems are related at short distances 
through the boundary condition corresponding to the 
short distance connection Si = S 2 



u' 2 {r c ) 
u 2 (r c ) 



1 R i 
l0g R~ 2 



Ri 



ui(r c ) 



(103) 



If we have only fixed the scattering length, the above 
condition becomes energy independent when the cut-off is 
small enough, which means that it can be evaluated with 
the zero-energy wave functions of the two-body systems 
1 and 2. Using the superposition principle, the previous 
zero energy wave functions can be written as 



ui(r) = vi(r) wi(r) , 

Oil 

u 2 (r) = v 2 (r) w 2 (r) . 

a 2 



(104) 
(105) 



These wave functions can be included in Eq. (|103|) , yield- 
ing the following relation between the scattering lengths 
ct\ and ol 2 of the two different problems 

A B D 
— = — + C + . 

ot\ a 2 Qi a 2 



(106) 



Therefore, if we make the hypothesis of Charge indepen- 
dence at short distances 4 



S n 



Sr, 



c 



(107) 



J np — "nn — •-'pp — "pp ■ 

and by making use of the superposition principle, we can 
write 

1 



uo, np {r) 

U ,nn(r) 

u o,p P (r) 



V0,np(r) 
v 0,pp( r ) - 



Oi0,np 
1 

1 



w ,n P (r) , 
■Wo,nn(r) ■ 



0,pp 



(r) 



0,pp 



(r)- 



1 



wo, P p(r) ■ 



C W 0,pp 
0,pp 



(r)-. 



(108) 
(109) 
(110) 
(111) 



Generally one might expect 



(1) 



S pP — S PP + a S P p 

. Our results are consistent with the expected smallness of the 
corrections. 
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NN 


LEP 


Renor. 


Reg. -I 


Reg. -II 


CD-Bonn [11] 


Exp. [11] 






(g^NN ~ o) 


{g^NN ~ 8) 


(g^NN ~ 20) 






np 


a [fm] 


input 


-23.737 


-23.738 


-23.738 


-23.74(2) 




r [fm] 


2.672 


2.678 


2.677 


2.671 


2.77(5) 


PP 


ao [fm] 


-17.806 


-18.350 


-20.088 


-17.46 


- 




r [fm] 


2.802 


2.799 


2.768 


2.845 


- 


PP(c) 


ao [fm] 


-7.706 


-7.824 


-8.265 


-7.8154 


-7.8149(29) 




r [fm] 


2.747 


2.641 


2.693 


2.773 


2.769(14) 


nn 


a [fm] 


-19.626 


-19.486 


-20.493 


-18.968 


-18.9(4) 




r [fm] 


2.771 


2.780 


2.763 


2.819 


2.75(11) 



TABLE II: NN low-energy parameters in the different scenarios. Renormalization only needs the np scattering length as an 
input parameter, all the other are calculated without ambiguities. The OBEP parameters have been fitted in the np case and 
kept the same in the other cases. Here (c) means Coulomb interaction is switch on. 





n+o-exch — 




Nijmegen >-« 




CD Bonn ■ ' 


np channel 




- m a = 490(17)MeV 
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nn channel 




a = -19.626 
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- pp channel 




' % = -7.706 
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FIG. 5: Renormalized phase shifts for the OBE potential with CSB OPE + a as a function of the cm. momentum in the 
singlet 1 S'o channel. In the left panel we show the fitted np phase-shift to the Nijmegen results Q. In the middle and right 
panels the predicted, pp(c) and nn are depicted and compared to the CD-Bonn result [111 ]. 



so we get the bilinear relations between all scattering 
lengths 



Br, 



A 



i'i> 



OLpp 

A c 
pp 



VP 



B,, 



B c 
pp 



D 



•mi 



c 



Dpp 
Olpp OLnp 

D c 

pp 

^pp &np 



(112) 
(113) 
(114) 



etc. In Fig. [4] we show the dependence of the scattering 
lengths as obtained from the np scattering length and the 
previous correlations. As can be seen, the correlations 
work rather well, confirming the idea that finiteness is 
a good criterion to implement charge independence of 
short distance operators. Numerical values are listed in 
Table [TTI for the experimental value of ao tnp . 

It is interesting to see how the short distance connec- 
tion works at finite energy and, in particular, if a given 
specific NN channel is able to predict the phase shifts for 



the remaining channels. In Fig. [S] we plot the extracted 
nn and pp(c) phase shifts when the OBE parameters have 
been fixed from the 1 >S'o Nijmegen np phase shifts. We 
have computed these phase shifts renormalizing in the np 
channel, i.e., fixing ao,np as input and integrating inward 
the Schrodinger equation, and then using Eq. (I107P we 
connect with the other channels. As we can see the short 
distance connection can be used to predict the 1 So phase 
shits for the rest of the channel with a high degree of 
accuracy. 



VI. CONCLUSIONS 

In this paper we have analyzed the charge dependence 
and charge symmetry breaking of the NN interaction. 
We have used the OBE model with exchange of 7r, a, to 
and p mesons and we have implemented CSB by means 
of pion mass splitting in the OPE potential and differ- 
ent nucleon masses. In particular, and as in previous 
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works, we have selected the 1 So np channel to fit scalar 
meson parameters, m a and ga-NN, as well as vector me- 
son couplings, g^NN and J p nn to the Nijmegen phase 
shifts [7|. A fine tuning problem arises when we using 
the customary regular boundary condition at the origin 
u(0) — 0. This problem appears in all np,nn,pp and 
pp(c) channels and large, ~ 40% , violations of SU(3) 
values of the g^NN coupling constant are needed. Tra- 
ditionally a great amount of effects such as multi-meson 
exchanges have been essential to explain the differences 
in phase shifts and threshold parameters for all np,n n,pp 
and pp(c) channels |3-fll| or the role played by p — uj |35l - 
[37j and/or ir — n mixing were invoked. These stan- 
dard approaches need very precise information on the 
interaction at all distances. 

However, once we admit incomplete knowledge of the 
interaction at short distances, it is possible to sidestep 
the problem of fine tuning by imposing a renormaliza- 
tion condition; at any stage of the calculation the scat- 
tering length is always kept fixed. This renormalization 
approach embodies short distance insensitivity. As a con- 
sequence, in the Charge Independent case, one can con- 
fortably take the experimental and/or SU(3) values for 
vector meson couplings. For the same reason we can only 
hope to quantitatively describe the relative changes due 
to the Charge Symmetry Breaking of the interaction at 
long distances. These considerations alone allow to ex- 
tract some universal information on the symmetry break- 
ing pattern where the np, nn and pp channels look very 
much the same at all energies even though the potentials 
are different and are indeed CSB. We have used a short 
distance condition to relate the renormalized np channel 
with the others nn, pp and pp(c). This short distance 
connection is so far an assumption based on finiteness 
but we have seen that reasonable results are obtained for 
low energy parameters and phase shifts. Our predictions 
for (AatciB, Arc/s) and (Aa C SB, Ar C ss) are compat- 
ible with the empirical one within the error estimation. 
This is in fact a remarkable result: all channels are gen- 
erated with just one scattering length, say np, and the 
long distance components of the potential where the CSB 
is, via physical pion and nucleon masses, explicitly built 
in. Our result is compatible with the interpretation that 
(relative) CSB sits at large distances. The absolute CSB 
is in a sense as uncertain as the short distance compo- 
nents of the NN force itself and cannot be determined 
independently of the Charge Invariant interaction. 
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Appendix A: Proton-proton fusion 

We would like to analyze further consequences of the 
short distance connection assumed by Eq. (I107p . An 
interesting process is the proton-proton fusion reaction 
pp — > d e + v e which is of central importance to stellar 
physics and neutrino astro-physics. In fact, it is the dom- 
inant solar neutrino source. The temperature in the Sun 
core is around T c = 15 x 10 6 K which means that we 
have protons of momentum p ~ (2m p T c ) 1 / 2 ~ 1.1 MeV. 
At these low energies, the reaction is dominated by the 
1 5'o — > d nuclear transition. The Gamow- Teller (GT) 
matrix element for this process (without MECs) is given 
by, 



A s M G t 



f 

Jo 



dr u 7 (r)M ,p P (r-) 



(Al) 



where uq jPP is the zero energy reduced wave function for 
the pp(c) system which can be related with the np prob- 
lem by Eq. (|T07l) . Then taking 

cto.np as input and inte- 
grating in we can calculate uq. pp . For deuteron we take 
as a first approximation the normalized bound state, 



i 7 (r) -> A s e" 



(A2) 



with 7^ = 0.2316fm~ and integrate inward the 
Schrodinger equation with negative energy E = 
—jj/M np . We obtain a value M GT = 5.189 fm to be 
compared to a more sophisticated one [34| using Argonne 
1^18 wave functions Mqt\avw — 4.859 fm. 
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In Fig. [6] we show the GT matrix element correlation 
with the np scattering length compared with the AV18 
calculation. Of course we have not included the tensor 
force which mixed S and D waves in the calculation of 



the deuteron. However we can appreciate that our num- 
bers are not very far from much more elaborate calcula- 
tions mi. 
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